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We present multiwavelength spectral analyses of two Fermi-LAT blazars, OJ 287 and 3C 279, 
that are part of the Boston University multiwaveband polarization program. The data have been 
compiled from observations with Fermi, Swift, RXTE, the VLBA, and various ground-based optical 
and radio telescopes. We simulate the dynamic spectral energy distributions (SEDs) within the 
framework of a multi-slice, time-dependent leptonic jet model for blazars, with radiation feedback, in 
the internal shock scenario. We use the physical jet parameters obtained from the VLBA monitoring 
to guide our modeling efforts. We discuss the role of intrinsic parameters and the interplay between 
synchrotron and inverse Compton radiation processes responsible for producing the resultant SEDs. 



I. INTRODUCTION 

Blazar jets are highly violent in nature and are 
dominated by ultrarelativistic particles. The SED 
of blazars consists of two spectral bumps. The low- 
energy component is due to synchrotron radiation 
emanating from relativistic particles, and the high- 
energy component (for leptonic jet model) is a result 
of Compton upscattering of the seed photon field by 
ultrarelativistic particles. The seed photons could ei- 
ther be the synchrotron photons produced in the jet 
(synchrotron self Compton, SSC) [tHHl, and/or ex- 
ternal disk photons entering the jet directly (external 
Compton disk, ECD) [H, Hj], and/or the disk photons 
getting reprocessed in the broad line region (BLR) 
(external Compton cloud, ECC) pi [ll|, and/or the 
dusty torus and then entering the jet. The spec- 
tral variability patterns and SEDs are important tools 
used for understanding the acceleration of particles 
and the time-dependent interplay of various radiation 
mechanisms responsible for the observed emission. 

Here, we analyze the multiwaveband SED of two 
Fermi-LAT blazars, OJ287 and 3C 279, using the 1- 
D multi-slice time-dependent leptonic jet model, with 
radiation feedback scheme of Q to gain understanding 
of the role of various intrinsic parameters and radia- 



tion processes in producing the resultant SEDs. 

We briefly describe the model of Q in fJTX| We dis- 
cuss our first results from this study in i jllll We sum- 
marize our results and give a brief description of future 
work in StTVj 



II. INTERNAL SHOCK MODEL 

The mode of acceleration of plasma electrons (and 
positrons) to highly relativistic energies and its loca- 
tion in the jet is still not completely understood. One 
way to comprehend the physics of particle accelera- 
tion is the internal shock model, in which the central 
engine (black hole + accretion disk) spews out shells 
of plasma with different velocity, mass, and energy. 
The collision between such shells gives rise to internal 
shocks (reverse (RS) and forward (FS)), which convert 
the ordered bulk kinetic energy of the plasma into the 
magnetic field energy and random kinetic energy of 
the particles. The highly accelerated particles then 
radiate and produce the emission observed from the 
jet. 

The collision of two plasma shells results in an emis- 
sion region as shown in Figure [TJ The treatment of 
shell collision and shock propagation is hydrodynamic 
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FIG. 1: Schematic of the emission region with RS traveling 
into the inner shell of bulk Lorentz factor (BLF) Ti, and 
FS moving into the outer shell of BLF with r o , such that 
Ti > r o . The primed quantities refer to the comoving 
frame and unprimed refer to the lab (AGN) frame. The 
comoving pressures p' fs and p' TS of the shocked fluids across 
the contact discontinuity (CD) are equal. A ra and Af s are 
the widths of the inner and outer shell after the collision 
obtained from the shock dynamics [l2| ]. 



and relativistic in nature [12j . 

The evolution of electron and photon population 
inside the emission region are governed, respectively, 
by, 



FIG. 2: Schematic of the radiative transfer in between 
the slices using the appropriate photon escape probability 
function. The unprimed n p h,~f w d/back values represent the 
photon densities in the forward and backward direction, 
respectively, in the comoving frame of the emission region. 
The rest of the unprimed quantities refer to the lab frame. 



which it gradually becomes transparent to radio fre- 
quencies is not simulated, as that would require the 
introduction of several additional, poorly constrained 
parameters. 



Multi-slice Radiation Transfer Scheme 



dn e (j,t) 
dt 

and 

dn ph (e,t) 
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A cylindrical emission region is considered to calcu- 
late the resultant spectrum in a time-dependent man- 
n e.\1i ner. The inhomogeneity in the photon and particle 
density throughout the emission region is realized by 
dividing the region into multiple slices, as shown in 
Figure [2 

The photon density of a zone and the probability of 
escape for a photon from that zone in a particular di- 
rection (forward, backward, or sideways), P, are used 
to calculate photon escape rates in that direction, ac- 
cording to 



Here, (d'j / dt)\ OBS is the radiative energy loss rate, 
due to synchrotron and SSC losses, for the electrons. 
Q e (7,i) is the sum of external injection and intrinsic 
7 — 7 pair production rate and i c ,csc is the electron 
escape time scale. n. p h, C m(e, t) and n P h,abs(e,i) are 
the photon emission and absorption rates correspond- 
ing to the electrons' radiative losses, and i p h,esc = 
(3/4)Rb/c is the photon escape timescale. The evolu- 
tion of the electron and photon population is followed 
in a time-dependent manner inside the emission re- 
gion and radiative energy loss rates as well as photon 
emissivities are calculated using the time-dependent 
radiation transfer code of [7|. 

The model follows the evolution of the emission 
region out to sub-pc scales and simulates only the 
early phase of 7-ray production. During this time, 
the radiative cooling is strongly dominant over adia- 
batic cooling and the emission region is highly opti- 
cally thick out to GHz radio frequencies. Thus, the 
simulated radio flux is well below that of the actual 
radio data. Also, the phase of the emission region in 



dn ph 

,fwd/back/sidc 



dt 



t 



ph,csc 



fwd/back/sidc 



(3) 

where i p h,csc is the photon escape timescale for a cylin- 
drical region. 

We use the scheme presented in equation [3] to cal- 
culate the radiation transfer within each slice and in 
between the slices Q. 



III. FIRST RESULTS 

The time-dependent model, calculating synchrotron 
and SSC radiation processes, developed in [7| has been 
used to reproduce the observed SED of two Fermi- 
LAT blazars, OJ287 and 3C 279 that are part of 
the Boston University multiwaveband monitoring pro- 
gram. We have collected the data from observations 
with Fermi, Swift, RXTE, the VLBA, and various 
ground-based optical and radio telescopes to construct 
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the multiwaveband SEDs of OJ287, and 3C 279, the 
latter of which corresponds to the quasar's optical 
high state as observed on 15 January 2006. 

OJ287 is a BL Lac object purported to have a black 
hole binary system [lj|. The blazar has exhibited 
spectral and polarization variability in the past [13]. 
Here, we apply the model of 0] with the ECD com- 
ponent included (Joshi et al., 2012, in prep.) with 50 
slices in the forward and 50 slices in the reverse emis- 
sion regions to analyse the SED of OJ287 as observed 
on 10/28/2008. Figure [3] shows the instantaneous and 
time-integrated simulated SED of OJ287 for that day. 
The model independent parameters [3] estimated from 
the SED, VLB A observations and variability on 
1-day timescale were used to develop an initial set of 
input parameters: 



5 « 16.5 

R « 3.5 x 10 15 cm 

B « 0.52 4 /7 G 

7 min « 1.1 x 10 3 

7max « 3.0 X 10 4 

q « 4.2 

6 ohs a 3.2° (4) 

Here, 5 is the Doppler boosting factor assumed to 
be equal to the BLF of the emission region (as esti- 
mated from VLB A observations). The symbols R & 
B represent the values of comoving radius and mag- 
netic field in the emission region and es, which is not 
a part of the model of 0] , is the ratio of magnetic field 
and electron energy density assumed to be equal to 1 
here. The 7 m i n & 7 max refer to the low and high en- 
ergy cutoffs of the electron energy distribution. The 
spectral index of the electron population is given by 
q and 6> b s is the observing angle inferred from VLBA 
observations. 

The initial set of parameters was modified to repro- 
duce the state of OJ287 as observed on 10/28/2008. 
Table |T] lists the parameters used for obtaining the re- 
sultant SED of the source shown in Fig. [3] These 
parameters result in a T s h ~ 16.4 for the entire 
emission region, B « 5.0 G and 7 max ~ 1.5 X 10 5 
for both forward and reverse emission regions, and 
7min, fs ~ 5.6 x 10 2 & 7 m i ni rs « 1.0 x 10 3 for forward 
and reverse emission regions, respectively. 

As can be seen from the figure, the lower-energy 
bump of the time-integrated simulated SED passes 
very close to the IR, optical, and UV data points, indi- 
cating that the synchrotron component is responsible 
for this part of the jet emission. The spectral upturn 
takes place in the soft X-rays at > 0.14 keV due to the 
presence of the SSC component in the simulation and 
the lower-energy part of the SSC component passes 
close to the X-ray data implying the dominance of 
SSC component in producing this part of the high- 
energy bump. The model (SSC+ECD), at this point, 



undcrpredicts the 7-ray photon flux suggesting that 
the contribution from the BLR might play a domi- 
nant role in reproducing this emission. 




FIG. 3: Simulated instantaneous and time-integrated (av- 
eraged over 1 day) SED of OJ287 for 10/28/2008. 

The flat-spectrum radio quasar (FSRQ), 3C 279, 
was observed in its optical high state on 01/15/2006, 
almost a month before it was observed by MAGIC 
emit ting in the TeV energy regime for the first time 
ever Figure [4] shows the instantaneous and time- 
integrated simulated SED of 3C 279 for that day. The 
model independent parameters 0] that were estimated 
using the SED, VLBA observations Q, and variability 
on 1-day timescale were used to develop an initial set 
of input parameters: 

5 w 15.5 

R w 2.5 x 10 16 cm 

B w 0.82 4 /7 G 

7min « 9.0 x 10 2 

7max « 2.2 x 10 4 

q w 4.3 

Oobs » 2.1° (5) 

All symbols refer to the same quantites as explained 
above and the entire emission region, as mentioned 
above, has been divided into 100 slices to analyze 
the observed SED of the source. The initial set of 
parameters was modified to reproduce the state of 
3C 279 as observed on 01/15/2006. Table U lists the 
parameters used for obtaining the resultant SED of 
the source shown in Fig. |4j These parameters re- 
sult in r s h ~ 16.6 for the entire emission region, 
B w 4.0 G and 

7max ~ 2.2 x 10 5 for both forward 
and reverse emission regions, and 7 m i ni f s ~ 5.8 x 10 2 
& 7min. rs ~ 1.3 x 10 3 for forward and reverse emission 
regions, respectively. 

As can be seen from the figure, the observed SED for 
01/15/2006 shows a high-energy bump that is indica- 
tive of a dominant SSC component and a suppressed 
EC component. The time-integrated simulated SED 
passes very close to the IR and optical data points, in- 
dicating that the synchrotron component responsible 
for the lower energy bump of the SED. The spectral 
upturn takes place in the soft X-rays at > 0.11 keV 
due to the presence of the SSC component in the sim- 
ulation. The lower-energy part of the SSC component 



eConf C110509 



4 



201 1 Fermi Symposium, Roma., May. 9-12 



reproduces the X-ray data quite well, suggesting the 
dominance of SSC component in producing this part 
of the high-energy bump. 




FIG. 4: Simulated instantaneous and time-integrated (av- 
eraged over 1 day) SED of 3C 279 for 01/15/2006. 



TABLE I: Model parameters used to reproduce the state 
of OJ287 and 3C 279 as observed on 10/28/2008 & 
01/15/2006, respectively. 



Source 


L w 


r, : 


r D 


q 




SB 






[10 47 ergs/s] 








[io- 2 ] 


[IO" 3 ] 


[IO" 2 ] 


3C279 


5 


35 


10 


4 


9 


3 


2.5 


OJ287 


4 


25 


12 


4.2 


9 


14 


1.0 


Source 


R 


#obs 


idisk 


Mbh 


Tja.cc 




[10 16 cm] 


[deg] 


[10 44 ergs/s] 


[10 8 Af S oi] 


[IO" 2 ] 


3C279 


3.7 


2.5 








OJ287 


4.0 


2.5 


2 


2 


6 



L w : luminosity of the injected electron population in the 

blob, I\ : BLFs of the inner and outer shells before 
collision, q: particle spectral index, e e : ratio of electron 
and shock energy density, eg: ratio of magnetic field and 

shock energy density, £ e : fraction of accelerated 
electrons, R: comoving radius, 8 bs'- viewing angle, Ldisk: 
accretion disk luminosity, Mbh: Mass of the BH, and 
7/ acc : accretion efficiency 



IV. DISCUSSION AND FUTURE WORK 

The time-integrated SEDs of both OJ287 & 3C 279 
need further adjustments of parameters that are listed 
in Table Q] in order to obtain a satisfactory fit. The 
external Compton component due to photons enter- 
ing the jet from the BLR and dusty torus needs to be 
incorporated in the existing model of Q to correctly 
reproduce the SEDs of blazars, especially for flat spec- 
trum radio quasars (Joshi et al. 2012, in prep.). 

Further, we plan to incorporate the effects of mag- 
netic field orientation, as inferred from polarization 
monitoring programs, on the resultant spectral vari- 
ability and SEDs of blazars. This would further aid us 
in the study of intrinsic parameter differences between 
various blazar subclasses, arising from the orientation 
of the magnetic field in the jet. 

We plan to study the evolution of SED of blazars 
from quiescent to flaring state in the light of the mod- 
ified time-dependent model of [Tfl (Joshi et al. 2012, in 
prep.). We will then compare the results with multi- 
waveband data gathered by the Fermi-LAT and other 
telescopes. 
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